Quantum phase transitions across a p-wave Feshbach resonance 
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We study a single-species polarized Fermi gas tuned across a narrow p-wave Feshbach resonance. 
We show that in the course of a BEC-BCS crossover the system can undergo a magnetic field- 
tuned quantum phase transition from a p^-wave to a p x + ip y -wa,ve superfluid. The latter state, 
that spontaneously breaks time-reversal symmetry, furthermore undergoes a topological p x + ip y to 
Px + ip y transition at zero chemical potential p. In two-dimensions, for fj, > it is characterized 
by a Pfaffian ground state exhibiting topological order and non-Abelian excitations familiar from 
fractional quantum Hall systems. 

PACS numbers: 03.75.Ss, 67.90.+Z, 74.20.Rp 
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Recent progress in controlling atomic interactions via 
an s-wave Feshbach resonance (FR) has led to a real- 
ization of an s-wave paired fermionic superfluid (SF) in 
a variety of degenerate atomic gases [1-3]. By tuning 
through FR such paired superfluid has been studied along 
the crossover between the BCS regime of weakly-paired, 
strongly overlapping Cooper pairs, and the BEC regime 
of tightly bound, weakly-interacting diatomic molecules. 

Although the study of this crossover has received con- 
siderable attention[4, 5], its equilibrium properties are 
already qualitatively well described by early seminal 
works[6, 7]. Furthermore, for a narrow FR, the thermo- 
dynamics can even be accurately computed analytically 
across the full range of the BCS-BEC crossover, with cor- 
rections controlled by the ratio of the FR width to Fermi 
energy [8]. 

In contrast, superfluids paired at a finite angular- 
momentum are characterized by richer order parame- 
ters (as exemplified by a p-wave paired He, heavy- 
fermion compounds, and d-wave high-T c superconduc- 
tors), and therefore admit phase transitions within the 
SF phase. Recent observation of a p-wave FR[9, 10] in 
i0 K and 6 Li, prepared in identical spin states, demon- 
strates a promising mechanism for a realization of p-wave 
superfluidity, attracting recent theoretical attention [11- 
13]. Because Pauli exclusion principle prevents identical 
fermionic atoms from scattering via s-wave channel, such 
p-wave scattering dominates. The tuning of the atomic 
interaction through such FR should allow access to afore- 
mentioned phase transitions and provides strong motiva- 
tion for our study. 

In this Letter, we study a p-wave superfluidity in a 
single-species polarized Fermi gas tuned across a p-wave 
FR. We show analytically that at low T, in the course of 
a BEC-BCS crossover this system can exhibit a second- 
order quantum phase transition as a function of FR de- 
tuning ujq (energy of the quasi-bound molecular state) 
from a p^-wave to a time-reversal breaking p x + ip y - 
wave superfluid. In three dimensions (3D), for a positive 
chemical potential /i > 0, the p x + ip^-wave superfluid 
exhibits gapless single particle excitations with momen- 



tum p = ±pj? z, that become gapped for \i < 0. Con- 
sequently, on quite general grounds we predict that in 
addition to above p x to p x + ip y phase transition, the lat- 
ter SF must furthermore undergo a topological quantum 
phase transition at loq(p = 0), from a gapped to gapless 
Px + ipy SF state[14] (Fig. l(a,b)) Associated thermo- 
dynamic changes (e.g., a heat capacity changing from 
activated to power-law form) should be observable. 
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1: (a) A schematic temperature-detuning phase diagram 



near a p-wave FR, for intermediate splitting of m = 0, ±1 FR 
(see below) . It shows a 2nd-order p x to p x + ip y SF transition 
at ujoc as well as a topological phase transition at o>o(m = 0) 
between two different p x +ip y phases, (b) , (c) show alternative 
phase diagrams for smaller and larger values of FR splitting, 
respectively. 



Moreover, when confined to 2D, for /i > 0, the p x + ip y 
superfluid is described by a fully-gapped Pfaffian (Moore- 
Read) ground state [15], that exhibits a subtle topolog- 
ical order characterized by excitations with non-Abelian 
statistics, familiar from fractional quantum Hall systems. 
At uj ([i = 0), this fully-gapped superfluid has been 
demonstrated to undergo a topological quantum phase 
transition to another strongly-coupled gapped p x + ip y - 
wave superfluid [16, 17]. 



2 



The accuracy of our theory is controlled by the ra- 
tio of the width of the FR (molecular decay rate) to ep. 
Because the low-energy decay rate of a finite angular mo- 
mentum diatomic molecule is suppressed by the centrifu- 
gal barrier, our theory can be made arbitrarily accurate 
by reducing the atom density, n. 

We start with a Hamiltonian of spinless fermionic 
atoms, resonantly interacting through a diatomic molecu- 
lar state with internal angular momentum ( "spin" ) I = 1: 



+ E 7f p* a p+§ a - P+ f + H • 



(i) 



Here ap is the creation operator of an atom with mo- 
mentum p and kinetic energy e p = p 2 /2m. fet cre- 
ates a molecule with momentum q, kinetic energy e q /2, 
and a labels its spin polarization. V is the volume of 
the system. e a controls the molecular rest energy re- 
lated to FR detuning, that can in general be polarization 
dependent [9]. The fermion- molecule matrix element g p 
is a constant go for p«A, quickly dropping off to zero 
for p ^> A. The ultraviolet (UV) cutoff A is associated 
with the smallest length scale in the problem, the size of 
the molecular (closed-channel) bound state. 

We first study elastic resonant scattering of two atoms 
with momenta p and — p into states p' and — p'. This 
will enable us to relate the bare parameters e a appearing 
in H to the experimental FR detuning u> a . In vacuum 
this process can only proceed through an intermediate 
molecular bound state, created by b a0 , at zero momen- 
tum and energy E = p 2 /m. The scattering amplitude 
/(p,p') is determined by the 4-fermion amputated cor- 
relation function. It can be deduced from the molec- 
ular Green's function at zero momentum, D a p(E) = 
8 a fj [E — e a — H(E)] 1 , expressed in terms of the molec- 
ular self-energy (fermionic polarization "bubble") 11(E), 

given by 11(E) = J ^A- f J^'S P r- 

Evaluating the integrals (with 3> E), we arrive at 
an exact expression for the p-wave scattering amplitude 



/(p,p') = E 



-(mg 2 /2ir)p aPa 



%(l + c 2 )-e a + Cl +i^p* 



(2) 



with ci = j±? J °° dk k 2 g 2 and c 2 = ^ J °° dk g\ quasi- 
molecular size (UV cutoff) dependent constants. 

Firstly, we observe that for a rotationally-invariant 
case of e a = e 0} Eq. (2) reduces to the I = 1 (p- 
wave) case of the most general expression for the scat- 
tering amplitude, fe(p) = [p~ 2l Fi>(p 2 ) — ip] _1 , for scat- 
tering from a central potential in angular momentum 
channel £, obtained simply from unitarity and analyt- 
icity. Fi(p 2 ) is a real function Taylor-expandable in pow- 
ers of its argument [18], that in our case can be read 



off from Eq. (2) and given exactly (for e\ 3> E) by 
Fx(p 2 ) = -A + Bp 2 , where A = - (e - ci) is the 

inverse scattering volume, and B = — ^ g -i (1 + C2). 

Hence, H in Eq. (1) correctly reproduces the low- 
energy two-body physics, regardless of the details of real 
interatomic interactions. The physical detuning position 
of the FR (corresponding to the molecular state that 
is real/virtual when it is negative/positive) is given by 
the low-energy poles of the scattering amplitude /(p, p'), 
■- , expressed in terms of the bare Hamiltonian 



- 1 + C 2 

parameters e a . In contrast to the s-wave case, where 
such renormalization only shifts the position of the FR, 
here it also leads to a multiplicative renormalization. 

We now proceed to the main topic of the Letter, 
namely to the study of thermodynamics of a fermionic 
atomic gas at a finite total density n, near a p-wave FR 
described by H, Eq. (1). It is convenient to work in the 
grand-canonical ensemble, = H — fin , with the chem- 
ical potential n determined by the average total atom 
density n = ± £ P (j2 a 2 ( b ap b ap ) + ( fit a P )) ■ 

cannot be diagonalized exactly. However, as we will 
show below, at T = the system is in a paired SF state, 
characterized by a finite ( b a Q ) = y/V B a for any value of 
detuning, uj a . This can be easily seen by noting that for 
a large positive detuning uj a 3> 2ep, there is a large gap 
to creating (quasi-) molecules, b a . Therefore they can be 
integrated out to lead to an effective BCS model, with 
an attractive p-wave four-Fermi interaction of strength 
Xbcs = Qo/i^a — 2fi), with virtual molecules mediating 
the attractive interactions between fermions. The result- 
ing finite Cooper-pair field ( a_ p a p ) in this BCS regime 
then always induces a finite (albeit exponentially small) 
(b a o) through a resonant go coupling. In the opposite 
BEC regime, io a < 2cf, we will show that \i decreases, 
tracking closely o> Q /2, and the energy can be lowered by 
binding a fraction of atoms into Bose-condensed p-wave 
molecules, stabilized by Pauli-blocking. 

Hence, to make progress we approximate the molecular 
field b ap by its condensate value ( b a o ) = VV B a , neglect- 
ing bosonic excitations. This reduces H p to a quadratic, 
mean-field BCS form (with B a playing the role of the 
Hubbard-Stratonovich field, i.e., the gap function), that 
can be easily diagonalized by a Bogoliubov transforma- 
tion. 

The key point is, that, independent of uj a , fluctuations 
about this molecular saddle-point approximation are sup- 
pressed in powers of a dimensionless coupling 

V2 
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(3) 



proportional to the ratio of the FR width to Fermi en- 
ergy. 7 can always be made sufficiently small by simply 
reducing the density n = (2mei?) 3 ^ 2 /(67r 2 ) [19]. Hence, 
for 7 <C 1, total number of excited molecular bosons 
S Q q^o ( ^aq^aq ) i s sma ll, and above molecular mean- 
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field approximation is quantitatively accurate, with cor- 
rections 0(7). In the BCS regime, this reduces to the 
standard criterion for the validity of the BCS mean- 
field corresponding to the gap much smaller than ep, 
namely Cooper pairs that are large compared to inter- 
atomic spacing[20]. However, in contrast to the stan- 
dard BCS model, the 7 <C 1 criterion is far less stringent 
and allows a quantitatively accurate analysis even in the 
molecular BEC regime, where from the BCS prospec- 
tive the atoms are strongly coupled with /i much smaller 
than the atomic gap. In this regime the limitation of 
the molecular mean-field theory simply corresponds to a 
small depletion of the molecular condensate. 

Hence for 7 <C 1, we can accurately compute the free 
energy F(B a ). As usual, its minimization with respect 
to the complex order parameter B, dF[B]/dB^ = leads 
to a BEC-BCS gap equation, 



(e Q - 2/i) B a = B p J 



d 3 k ggfcqfeg / E k , 
(2tt)3 E k tanh \2k B T>' (4) 



whose solution, together with the total atom conservation 
equation 



n = 2n b + n a 



(5) 



\2k B T ) 



IS 



determines the molecular condensate amplitude, B and 
the chemical potential, /i as a function u a , T and n. 

In above, n a {pt) = \ J -0^ 1 - tanh 

the number of open-channel BCS-paired atoms, with 

E k = \J _ ^j 2 _)_ 4g|| B • k| 2 the Bogoliubov atomic 

quasi-particle spectrum, and rib = 5Z Qq (&a q k*q) the to- 
tal number of closed-channel molecules. 

The complex vector order parameter B encodes the 
three molecular states ip m with spin projection m = 0, ±1 
through the relations ipo = B z , -01 = — (B x + iB y ), 
= -^j (B x — iB y ) [18]. It can be parameterized in 
terms of two real vectors 



B = u + iv, 



(6) 



that by a gauge transformation B — > e IV B can always 
be conveniently chosen to be orthogonal, with u • v = 0. 
Consequently, if either one (but not both) of the vectors 
vanishes (e.g., v), then the projection of the angular mo- 
mentum on the axis parallel to the other vector (u) is 
equal to zero. If, on the other hand, u = v (with u = |u| 
and v — |v|), then the projection of the angular momen- 
tum on v x u is equal to 1 [18]. 

Focusing first on the T = case, by integrating Eq. (4) 
we find the free energy in this u, v representation 



^"'^ = V (u 2 a + v 2 a ) [w a -2 f i + a 1 In {a (u + v)}} 

(7) 



+a 



C2 

u a + V 3 



U + V 



(u 2 + v 2 y + ^(u 2 -v 2 y 



with 01 = (it^VmAf) M, «2 = (f Tfe7) £F/n, 

ao = e 5 / 6 (/i/ei?) 1 / 4 (7/8n) 1 / 2 . As usual, the logarithmic 
term arises from the Fermi surface (k 2 /2m ~ ji) con- 
tribution to the integral of Eq. (4). It is only a good 
estimate for [i > when the logarithm itself is large, 
and vanishes otherwise. On the other hand, the quar- 
tic terms a 2 [ (B* • B) 2 + | |B • B| 2 ] in F arise from mo- 
menta k <~ A. 

Minimization of F(u, v) leads to the phase diagram 
as a function of detuning oj a , displayed in Fig.l. For 
luq 3> 2tp (BCS regime), the molecular density is van- 
ishingly small and atom conservation, Eq. (5), constrains 
[i to be close, but slightly below ep. The minimum of 
F(u, v) is then controlled by the balance of the quadratic 
and logarithmic terms in u and v, with the quartic 
terms exponentially small. In the isotropic case with 
molecular spin-independent FR, uj a = ujq, F is invari- 
ant under simultaneous rotation of u and v, and ex- 
hibits three (plus symmetry related) local extrema: (i) 
normal state with u = v = 0, (ii) p^-wave supcrfluid 
with u ^ 0, v — 0, and (iii) p x + ip y -wave superfluid 
with u = v 7^ 0. The absolute minimum of F is given by 
(iii) with u = v w e~^°~ 2eF ^ ai , and the condensa- 
tion energy ratio TZ — F Px+ i Py /F Px w e/2 > 1, robustly 
giving p x + ip y as the ground state in the BCS regime, 
consistent with the A\ supcrfluid ground state of a fully 
spin-polarized 3 Hc.[ll, 21, 22]. 

In the complementary loq < 2ep crossover regime, it 
becomes favorable to pair-up a finite fraction of atoms 
into closed-channel Bose-condensed molecules at energy 
uoq. This forces \i to track wo/2 below tp. Consequently, 
the condensate density n& w |B| 2 is no longer exponen- 
tially small, growing to a finite fraction of the overall 
density n, and as a result the logarithmic term in Eq. (7) 
becomes negligible. The condensate density |i?| 2 is then 
primarily determined by the atom number Eq. (5), with 
the gap equation simply constraining \i to track 
to the accuracy 0(j). In the isotropic case, it is clear 
from Eq. (7) that although the quadratic B* • B and 
the first quartic (B* • B) 2 terms are invariant under ro- 
tation in the u — v space, the positive coefficient of the 
|B • B| 2 = (u 2 — v 2 ) 2 quartic term guarantees that F(u, v) 



is again minimized by u 



\y / n - n a (u> Q /2), cor- 



responding to the time-reversal-breaking p x + ip y -w&ve 
superfluid ground state. In this idealized isotropic case, 
the finite T transition from the normal to the p x + ip y 
supcrfluid state is described by the complex 0(3) vector 
model, F = tB* ■ B + Ai (B* • B) 2 + A 2 |B • B| 2 . 

However, in the experimental realization, the mag- 
netic dipole interaction generically splits the FR into 
two m = and m = ±1 resonances [9]. Choosing x-axis 
along the external magnetic field, this can be captured 
in our model Eq. (1) by taking uj x = ojq — 5 ^ u) y = 
lu z = wo- This orbital easy axis (experimentally S > 0) 
anisotropy substantially enriches the phase diagram aris- 
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ing from minimization of F[u, v], Eq.7. A technically 
involved analysis [8] predicts for an intermediate FR split- 
ting 1/(1 + c 2 ) < 5 1 {run) < 5/(4c 2 ) (Fig. 1(a)), that for 
luq < uj 0c the free-energy is minimized by a u = ux, 
v = state corresponding to a p^-wave m = super- 
fluid. At low T, for uj > ^Oc we instead find that the 
free-energy is minimized by a u ^ f / 0, characterizing 
a state that breaks time-reversal symmetry and is a co- 
herent superposition of p x and p x + ip y supcrfluids (that 
we refer to as simply a p x + ip y state). It is easy to see 
that the quantum phase transition from the p x to the 
Px + ipy superfluid is second-order. The global (gauge) 
U(l) symmetry is spontaneously broken in the p x super- 
fluid, and rotational 0(3) symmetry is explicitly broken 
(down to 0(2)) by the dipole interaction that acts like 
Ising (easy axis) anisotropy. Hence it can be shown [8] 
that the T = p x to p x + ip y transition at lo$ c is in the 4D 
0(2) (XY) universality class. The corresponding critical 
point controls the shape of the low T phase boundary 
with T c i(cjo) ~ |o>o — wo c | Z!y , and z = 1 and v = \. 

From the general structure of F near T c2 (<~ ep for 
luq < 2ep and exponentially small in the BCS regime), 
where the nonanalytic logarithmic term in Eq. (7) is sup- 
pressed (through the hyperbolic tangent in Eq. (4)), it is 
possible to show that for an arbitrary small FR splitting 
6 > the normal-superfluid transition is always into the 
p x state (with the transition in the 3D XY universality 
class), in contrast to the isotropic case, 6 = 0, where it is 
always into the p x + ip y superfluid [8]. 

For a sufficiently small FR splitting 5/(na 2 ) < 1/(1 + 
c 2 ) (Fig. 1(b)), the T = phase transition disappears 
and at low T the p x + ip y superfluid extends to all detun- 
ing, undergoing a transition to the p x state at T cl (w ) 
(~ ep for luq < 2ep) only near the normal-to-supcrfluid 
T C 2 (wo), below which the p x superfluid is confined in a 
sliver of T cl < T < T c2 set by In [T c2 /T cl ] ~ 8. This is 
consistent with the isotropic case 8 = 0, where p x + ip y 
state extends all the way up to T c2 . In contrast, for a 
large FR splitting (Fig. 1(c)), 8/(na 2 ) > 5/(4c 2 ), it is 
clear from F, Eq. (7) that the dipole anisotropy domi- 
nates all other terms and leads to the p x state over the full 
T — ivo phase diagram. However, it is important to note 
that for a very large FR splitting, 6, given long energy 
relaxation times, near the m = 1 resonance the atomic 
gas may remain in a metastable p x + ip y state, despite 
the existence of a lower energy p x state. 

The analysis presented in this Letter establishes the 
stability of p x + ip y SF phase over a large portion 
of the phase diagram, in an experimentally realizable 
system in which the chemical potential p can be di- 
rectly tuned through zero by adjusting FR detuning 
luo. This allows us to take advantage of the established 
literaturc[6, 14, 16, 17], that demonstrates that such SF 
state undergoes a topological phase transition at \i = 
between two types of p x + ip y SFs. In 3D these are 
distinguished by the existence of gapless atomic excita- 



tions for \x > (due to nodal points at p = ±Pf z) and 
a fully gapped spectrum for ji < 0. In 2D (to which 
our computation of F can be easily extended without 
any qualitative changes [8]) both of these states are fully 
gapped but nevertheless are distinguished by topological 
order[16, 17]. Thus our work explicitly shows that de- 
generate Fermi gases near a p-wave FR give a concrete 
realization of non-Abelian topological states, which are of 
interest in strongly correlated systems with possible ap- 
plications to fault-tolerant quantum computing [23]. De- 
tecting phases of p-wave superfluids and transitions be- 
tween them (illustrated in Fig. 1) in cold atomic gases 
remains an open problem. 

We acknowledge support by the NSF DMR-0321848 
(L.R.), DMR-9984002 (A.A.), and the Packard Founda- 
tion (A. A., L.R.), and thank J. Bohn, and D. Jin for 
discussions. 

Note added: After this paper was posted, a manuscript 
by Chi-Ho Cheng and S.-K. Yip [24] appeared that stud- 
ied the same problem. For large (5/{na 2 ) > 5/(4c 2 )) and 
small (8/{na 2 ) < 1/(1 + c 2 )) Feshbach resonance split- 
ting their results coincide with ours (summarized by our 
Figs. 1 (b,c)). However, results disagreed (see Fig. 1(a) in 
version 2 of our manuscript) for the narrow range of the 
intermediate values of the Feshbach resonance splitting, 
1/(1 + c 2 ) < 5/(na 2 ) < 5/(4c 2 ), for reasons at the time 
not understood by either groups. Upon further analy- 
sis we found that Cheng and Yip's predictions [24] for 
this narrow range of FR splitting was correct. Our mis- 
take stemmed from not appreciating the dominant (in 
the narrow intermediate range of FR splitting) contri- 
bution to the atom number coming from the long tail 
in the atomic momentum distribution function. Correc- 
tion of this mistake is reflected in the updated version of 
Fig. 1(a), but does not affect any other result in our origi- 
nal manuscript. We credit Cheng and Yip for the correct 
version of Fig. 1(a). Detailed analysis will be presented 
in the manuscript currently in preparation [25]. 

We regret that due to our oversight this note did not 
make it into the version of our manuscript published in 
Physical Review Letters. 
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